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ABSTRACT 

We study the density of clustered immobile reactants in the diffusion-controlled single 
species annihilation. An initial state in which these impurities occupy a subspace of 
codimension d leads to a substantial enhancement of their survival probability. The 
Smoluchowski rate theory suggests that the codimensionality plays a crucial role in 
determining the long time behavior. The system undergoes a transition at ci = 2. For 
< 2, a finite fraction of the impurities survive: riiit) = nj{oo) + const x log(t) / \/t for 
d = 2 and nj{t) = n/(oo) + const/-\/t for d > 2. Above this critical codimension, d > 2, 
the subspace decays indefinitely. At the critical codimension, inverse logarithmic decay 
occurs, ni{t) ~ log(t)~"'^'^''^''. Above the critical codimension, the decay is algebraic 
n/(t) ~ j^-"(<^''^). In general, the exponents governing the long time behavior depend on 
the dimension as well as the codimension. 
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I. Introduction 

The kinetics of diffusion-controlled chemical reactions have attracted much interest 
recently. For simple homogeneous reaction processes, substantial theoretical knowledge 
is available |]l], |, ^, H, ||, P]. For single species bimolecular reactions it is well known 
that for d <2, spatial correlations between reacting particles are important in the long 
time limit, while for d > 2 reactants are effectively transparent, and a universal decay of 
the density occurs. A number of reaction processes, such as the annihilation processes, 
A + A ^ and the aggregation process, A + A ^ A belong to this universality 
class. 

A recent generalization of this process to heterogeneous situations was shown to 
exhibit a rich array of asymptotic behavior. When the reactants have a polydisperse 
distribution of diffusion coefficients, non-universal decay kinetics may occur. Such a pro- 
cess is well suited for describing reactions that involve particles with different masses. 
The case in which a small number of particles move according to one diffusion coeffi- 
cient and the bulk according to another is especially interesting because of its simplicity. 
The survival probability of a single "impurity" particle immersed in a background of 
identical particles, nj{t), depends in a nonuniversal fashion on the diffusivity of the 
impurity. In one dimension, the results are especially intriguing. While for the ag- 
gregation process, a mapping to a random walk in three dimension enables an exact 
solution 0, for the annihilation process only an approximate theory is available 0. 
The survival probability of an immobile impurity particle is equivalent to the fraction 
of unvisited sites by annihilating random walkers [^]. Additionally, as the annihilation 
process is equivalent to the T = Ising model with Glauber dynamics [l^, nj{t) equals 
the fraction of "cold" spins, i. e., the number of spins that did not flip up to time t. 
Numerically, both time-dependent and finite size simulations p|, |TI], [l^ , as well as time 
series studies [|], [1^ suggest that the impurity survival probability follows a non-trivial 
decay nj{t) ~ t^"' , with aj = 0.37, while the corresponding theoretical exponent for 
the aggregation case is = 1. In addition, for trimolecular annihilation, A-l-A-l-A —* 
it was found that the impurity decays faster than a power-law but slower than an ex- 
ponential ni{t) ~ exp ( — log(t)^/^) p. Hence, in one-dimension, the impurity decay 
kinetics are highly sensitive to microscopic details of the reaction process. 

In this study, we present a generalization of the isolated impurity problem. We 
examine the collective effects of "trapped" or equivalently immobile impurities by con- 
sidering the initial condition in which the immobiles are clustered. For simplicity, the 
immobiles occupy a subspace of dimension dj embedded in a d-dimensional space. For 
convenience we introduce the codimensionality d = d — dj. Note that the case d = d 
corresponds to the isolated impurity problem. We are primarily interested in the rele- 
vance of the subspace, if any, to the asymptotic behavior of the survival probability. In 
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the long time limit, neighboring immobile reactants "shield" each other, and as a result, 
the immobile reactants decay significantly slower than the background. The subspace 
becomes more pronounced and consequently, the subspace interface acts as an absorbing 
boundary to the background. The geometry of the system reduces to a (i-dimensional 
one. By applying the Smoluchowski rate theory to the background in dimension d and 
to the impurity subspace in dimension d, theoretical prediction concerning the asymp- 
totic form of the survival probability are made possible. The subspace undergoes a 
survival-extinction transition at d = dc = 2. Below this critical codimension, a finite 
fraction of the immobiles survive, while at a higher codimension they decay forever. 
For d = 1, the approach to the final density is an algebraic one ni(t) — n/(oo) ~ 1/Vt 
for d > 2, with a logarithmic correction at d = 2, ni{t) — n/(oo) ~ \og{t)/^/t. At the 
critical codimension, d = 2, an unusual logarithmic decay occurs n/(t) ~ log(t)~"^'^''^\ 
while for J > 2 a dimension dependent decay ni(t) ~ occurs as well. 

II. The Isolated Impurity Problem 

In this section we review the rate equation theory and apply it to the isolated 
impurity problem. This approach truncates the infinite hierarchy of equations describing 
the density at the first order. Above the critical dimension, spatial fiuctuations are 
asymptotically irrelevant and thus, this theory is exact. However, this theory can be 
extended to arbitrary dimension, and is especially attractive due to its simplicity. 

In the lattice version of the homogeneous single species annihilation process, particles 
hop independently with rate Db to any one of their nearest-neighbor sites. An attempt 
to land on an occupied site results in the removal of both particles from the system. 
Hence, the density, ^^(t), obeys the following rate equation, dns/dt oc —ubb, where 
ubb is the density of pairs of neighboring particles. This approach leads to an infinite 
hierarchy of rate equations and is of limited practical use. Alternatively, one assumes 
that the annihilation rate is proportional to the fiux experienced by a particle, duB /dt (x 
—jriB fl^. This fiux can be evaluated by placing an absorbing particle in a diffusing 
background of density ub- Since j is proportional to n^, the density is described by 
the following rate equation, 

with k'^B being the reaction rate. Above the critical dimension dc = 2 the reaction 
rate is constant, and the above theory is exact. Otherwise, the reaction rate is time 
dependent since the target particle is surrounded by a depletion zone the size of the 
diffusion length y/Dt. In appendix I, we detail a heuristic derivation of the reaction 
rate in arbitrary dimension using the quasistatic approximation. In general, kfj, the 
reaction rate between particles of diffusivities Dj and Dj has the following long time 
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form, 

{A^^{Di + Dj)/t d=l- 

kij^lA,iDj + Dj)/log{iDi + Dj)t) d = 2- (2) 

[ Ad{Di + Dj) d> 2. 

Note that both kfj and Dj are rates and that is a dimensionless prefactor. By 
introducing a modified time variable 

m= fkU{t')dt', (3) 
Jo 

the rate equation simplifies, dns/dz = —n\. The asymptotic solution to this equation 
is ~ 1/2;. Evaluating the modified time variable, we arrive at the following 
asymptotic behaviors of the homogeneous single-species annihilation process. 



1/Vt d=l] 
nB{t)^ {\og{t)/t d = 2- (4) 
1/t d>2. 

Interestingly, these results are asymptotically exact ^, despite the assumptions in- 
volved with the rate theory. Moreover, the same results are obtained for the aggregation 
process A + A A as well. In lower dimensions, the reaction proceeds with a slower 
rate, since particles are effectively repelling each other. On the other hand, in high 
dimensions spatial correlations are practically irrelevant and the density decays faster. 
The critical dimension is characterized by typical logarithmic corrections. 

In the impurity problem, the survival probability of a single impurity with diffusivity 
Dj, ni{t), is investigated. Eq. (|I]) can be generalized to situations with more than one 
species 

fir} T 

-k^jHsnj, (5) 



dt 

with k'j^j given by Eq. (0). Again, it is useful to rewrite the rate equation in terms of 
the modified time variable z defined by Eq. (^), drij/dz = —{k'^j/k'^^)nBni. Note that 
in the long time limit the rate ratio approaches a constant that depends only on the 
diffusivity ratio. By substituting the asymptotic form of the background density 
one finds a purely algebraic dependence of the impurity density rij ~ z^bi/^bb for d ^ 2. 
For the case d = 2 the leading asymptotic correction to the rate ratio is important and 
a slightly more detailed calculation is needed. Using the aforementioned forms of z and 
kfj, we find the following behavior for the impurity density. 



nj{t) ~ <j ^-{Ds+Dj)/2Ds iog(i)7 d = 2; (6) 
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Hence, the diffusivity ratio, r = [Dj + Db)/2Db, governs the long time kinetics. The 
logarithmic correction for the case = 2 is characterized by 7 = r(l— logr) |T5|. For the 
case Dj = 0, the decay exponent equals 1/ VS = 0.353 when d = 1, and 1/2 when d > 2. 
While the latter value is exact, the former is only approximate Interestingly, this 
approximate value is quite close to the observed numerical value 0.37. However, as this 
approximation is uncontrolled, its accuracy can widely vary, and for the exactly soluble 
aggregation process, the discrepancy in the exponents is much larger. To summarize, 
both in the supercritical regime and in the subcritical regime, the decay of the impurities 
depends on the diffusivity ratio in a nongeneric fashion. 

III. Kinetics of Clustered Impurities 

The following questions arise naturally from the above theory: Can the presence 
of neighboring impurities increase the survival probability of an impurity? If yes, to 
what extent? To answer these questions we start with an initial configuration where 
impurities occupy a subspace of dimension dj. The codimension d can be conveniently 
defined as d = d — dj. For example, when d = 1, the impurities initially occupy a line 
in 2D, a plane in 3D, etc. Figure |I] illustrates the initial configuration for the case d = 1 
in two spatial dimensions. Note also that the single impurity problem corresponds to 
the special case d = d. To preserve the geometrical properties of the impurity subspace 
as the reaction evolves, we restrict our attention to the case Dj = 0. 

Let us consider a line of impurities in 2D. The presence of nearby static particles can 
only increase the survival probability of an impurity, and thus, the survival probability 
is bounded by rij ~ t"^/^, the corresponding result for the single impurity case for d > 2. 
However, the background decay ub ^ is stronger and consequently, the subspace 
becomes more and more pronounced in the long time limit. Hence, the impurity- 
background interface is equivalent to an absorbing boundary for the mobile reactants. 
This absorber is not a perfect one since not all sites are occupied with impurities. 
However, it is well known that in the long time limit a partial absorber is equivalent to 
a perfect one. A depletion layer of width y/Dt develops around the subspace and the 
background density profile is strongly suppressed in this depletion zone. As a result, the 
geometry of the system is drastically altered. A line in 2D reduces to a one- dimensional 
geometry. In other words, the codimension becomes the relevant dimension. This simple 
conclusion has a striking effect on the long time kinetics of the impurities. 

One possible way to tackle this problem is to describe this spatial inhomogeneity 
by the reaction-diffusion equation with proper boundary conditions accounting for the 
impurity interface. However, the leading term in the reaction-diffusion equation is the 
diffusion term and this approach is equivalent in the long time limit to the quasistatic 
approximation, or namely, the Smolechowsky rate theory. This theory, appealing in its 
simplicity, leads to new and interesting behaviors for the clustered impurity problem. 



4 



o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


• 


o 


o 



Figure 1: A line of immobile impurities (bullets) in a two dimensional background 
(circles) . 



As discussed above, the relevant dimension is the codimension. Thus, we apply the 
d-dimensional rate equation to the impurities, 

^ = -kBlUBUi, (7) 

with the reaction rates of Eq. (^). The equation describing the background is left 
unchanged, and one simply substitutes the density from Eq. (^. In contrast with 
the previous analysis, introduction of a modified time variable would not simplify the 
algebra, since the equations involve different dimensions and consequently, different 
intrinsic time scales. Instead, the impurity survival probability, ni{t)/ni{0), is obtained 
by an integration of the above equation, 

nj{t)/nj{0) = exp (^- |J dt'kij{t')nB{t')^ . (8) 

Since the case d = d reduces to the single impurity case (see Eq. (P)), we concentrate 
on the case d < d only. By substituting the proper values for the (i-dimensional back- 
ground density and the (i-dimensional reaction rate into Eq. we find the following 
asymptotic impurity densities. 



ni{oc) + const x log(t)/v^ d = 1 and d = 2: 

71/ ( cxo ) + const /i/t d = 1 and d > 2: 

log(t)-^d72A, d = 2andd> d 

r^d"/2A, d> 2 and d>d. 



(9) 



This rich behavior follows directly from the annihilation rate, i. e., the integrand in 
Eq. (H). If k%j{t)nB{t) decays faster than the integral remains finite in the infi- 
nite time limit, and a finite fraction of the impurities survive the annihilation process. 
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Otherwise, the subspace vanishes. For J < 2, this integrand decays faster than 1/t, 
while for d > 2, the integrand is dominated by 1/t. Consequently, the system exhibits 
a survival-extinction transition aX d = 2. Below this critical codimension, a fraction of 
the impurities survive, while they decay indefinitely at a higher codimension. 

Interestingly, the approach to the final density, n/(t) — n/(oo) oc log(t)/v^ for a line 
in 2D {d = 1) is reminiscent of the single impurity decay in 2D. These cases differ in 
that ni{oo) does not vanish for the clustered case. They also differ in their logarithmic 
correction. The critical case is characterized by inverse logarithmic decay since the 
annihilation rate is proportional to 1/ (Uog(i)) = dlog(log(i)) /dt. Both the critical 
case and the supercritical regime follow a power-law with the exponent, a{d,d) — 
A^/2Ad. Thus, both decays depend on the dimension as well as the codimension. Only 
in the extreme case, d — d, corresponding to the single impurity case, cancels 
out and a universal decay is expected, n/ ~ Generally, a detailed calculation 

for the prefactors is necessary in order to find the various decay exponents. Since 
the decay in the clustered case d < d is slower than in the isolated case d = d, we 
learn that the prefactor A^ is an increasing function of the spatial dimension d. This 
observation is consistent with the fact that the initial reaction rate is given by kjjlt — 
0) = 2{Di + Dj)zd, with Zd the number of neighboring sites in d-dimensions. Indeed, Zd 
is an increasing function of d. Above the critical dimension d — 2, the effective reaction 
rate depends weakly on time. Hence, the assumption Ad oc Zd leads to an approximate 
value for the decay exponent a{d,d) = z^/2zd. For a simple square lattice one has 
the Zd = 2d or a{d,d) = d/2d. Note also that this approximation improves as the 
dimension and the codimension increase. 

The above results can be easily generalized to arbitrary dimensions. Such a gener- 
alization is nontrivial only when the temporal nature of the reaction rate is dimension 
dependent, or namely, below the critical dimension. Using the reaction rates of Eq. (A3), 
we evaluate the immobile impurity densities, 

' t-rf/2^+'/' d <2 and d = d; 

n/(oo) + const x t-(^-<^)/2 d <2 and d < d < 2; 
n/(oo) + const x \og{t)t~^'^-^y'^ J < 2 and d = 2; 

ni{t) ~ { n/(oo) + const x r(2-(i)/2 J < 2 and d > 2; (10) 

iog(t)(i+iog2)/2 J = 2 and d = 2 

log{t)-^^/^^i J = 2 and d >(i; 

[ t-A^/2A^ J > 2 and ci > J. 

To summarize, the impurity subspace survives only when d < 2 and d > d. The 
behavior for d < 2 is influenced by the background density behavior ~ i"*^/^. As a 
result, the approach towards the limiting density is algebraic with a vanishing decay 
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exponent {d — d)/2 = dj/2, for the case d < d < 2. It will be interesting to see how 
well the above results compare with Renormalization Group studies in the vicinity of 
the critical codimension d = 2 — e. 

For completeness, we briefly discuss the early time behavior of the system. We focus 
on the case where all lattice sites are initially occupied, such that nsiO) = ni{0) = 1. 
Following the above discussion, the initial background-background reaction rate is given 
by k'^si^ = 0) = 2Dzd- From Eq. (|lD, the background density is found 

neit) = {l + 2Dzdt)~^ t ^ 0. (11) 

On the other hand, only interfacial sites contribute to annihilation of impurities, and as 
a result k%j = D{zd — -^d-d)- By substituting this rate and the early time background 
density into Eq. (|^), the impurity density is calculated in the early time regime, 

ni{t) ^ (1 + 2Dzdtf^'^'^^ t 0. (12) 

The above exponent, P{d,d), equals the reaction rate ratio P{d,d) = {zd — Zd-djf^^d- 
This exponent should not be regarded as an asymptotic one, since it is relevant only 
for a short time. In addition to the dimension dependence, the early time behavior 
depends on the lattice structure as well. As the reaction process evolves, such details 
become irrelevant, and the general asymptotic behavior is recovered. Note that there 
is no sign of a critical codimensionality in the early stages since the system is still 
(i-dimensional. After waiting a sufficiently long time, the geometry changes, and the 
codimension governs the kinetics. 

IV. Conclusions 

We have studied the kinetic behavior of simple subspaces of immobile reactants in 
the annihilation process A + A — > 0. Asymptotically, such an inhomogeneity leads to 
a change in the geometry of the system, and the codimension becomes the relevant 
parameter. The Smolechowsky theory shows that a transition from survival to extinc- 
tion takes place aX d = 2. For d < 2, and d > d, a finite fraction of the impurities 
survive, while for d > 2 the impurity subspace eventually vanishes. Furthermore, the 
asymptotic behavior of the impurity density depends on the dimension as well as the 
codimension. 

This study suggests that the rate theory can be extended to heterogeneous situations. 
Moreover, a time- dependent reaction equation is equivalent in the long time limit to the 
detailed reaction-diffusion equation. The success of this theory is remarkable especially 
considering its simplicity. It will be interesting to apply the same mechanism to more 
complicated processes, such as multispecies reactions. 
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Appendix A. The Quasistatic Approximation 

The reaction rate, kfj, can be found by evaluating the flux j experienced by an 
absorbing test particle of radius R due to a diffusing background of concentration 
Cq. Hence, the diffusion equation dc/dt = DVc is solved under the initial conditions 
c\t=o = Co and the absorbing boundary condition c{r)\r=R = 0. Although an exact so- 
lution can be obtained, we present a simpler approximate technique. Since the length 
scale governing the diffusion process is VDt, one assumes that the absorber does not 
affect the background density at distances larger than ^/Dt. Inside this depletion zone 
the Laplacian term dominates the spherically symmetric diffusion equation, 

Dr'-'^r''-^^ = 0, R<r<VDi. (Al) 
or or 

The quasistatic approximation imposes an additional time-dependent boundary condi- 



tion cl^^^qyi = Co [|16|. The concentration profile is readily obtained for R < r < y Dt 



Co ((r/i?)2-^ - l) / ({y/Di/Rf-'' - l) d<2- 
c(r, t) ~ <j Co \og{r/R)/\og{VDi/R) d = 2; (A2) 

Co (l - (r/i?)2-'^) / (l - {./Di/Rf-'') d>2. 

Above the critical dimension d^. = 2, c{r,t) approaches a limiting density profile. The 
reaction rate is given by the calculating the total flux seen by the test particle j = 
DSdR'^~^dc/dr, with Sd the surface area of the ci-dimensional unit sphere. We quote 
the leading asymptotic term of the reaction rate, k'^ = J/cq, 

k oc I D / log{Dt) d = 2] (A3) 
[ DR'^-^ d>2. 

In case where the target diffusivity Dt is nonzero, the effective diffusion constant is 
D + Dt- We derived the continuum rates, however, the lattice counterparts can be 
conveniently obtained by setting R=l. 
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